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We investigate the thermodynamics of a one-dimensional Hubbard model with bond-charge in- 
teraction X using the transfer matrix renormalization group method (TMRG). Numerical results 
for various quantities like spin and charge susceptibilities, particle densities, specific heat and ther- 
mal correlation lengths are presented and discussed. We compare our data also to results for the 
exactly solvable case X/t = 1 as well as to bosonisation results for weak coupling X/t <^ 1, which 
shows excellent agreement. We confirm the existence of a Tomonaga-Luttinger and a Luther-Emery 
liquid phase, in agreement with previous studies at zero temperature. Thermal singlet-pair correla- 
tion lengths are shown to dominate density and spin correlations for finite temperatures in certain 
parameter regimes. 
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I. INTRODUCTION 

The Coulomb interaction between electrons in a solid 
does not only lead to diagonal terms like the Hubbard 
interaction. It is also responsible for other non-diagonal 
terms that are often neglected since the corresponding 
coupling constants are assumed to be small. However, 
this is not always true. Hirsch (see e.g. refs. and 
references therein) has argued that especially the bond- 
charge interaction X may be responsible for interesting 
physics based on the breaking of particle-hole symmetry. 

In refs. it has been shown within a BCS-type 

analysis that the bond-charge interaction can induce an 
effective attraction between particles at densities n > 1, 
which can even lead to the formation of Cooper pairs 
and superconductivity. More evidence for this sce- 
nario has been obtained in analyticalSiLSiSiiSiiLiS and 
numericaliii^iiSii& studies of one-dimensional variants 
and an effective model for Cu02 planesii. Other ma- 
terials where bond-charge interaction might play an im- 
portant role are polyacetylene (see ref. and references 
therein) and the Bechgaard saltsiS. 

The Bechgaard salts show a rich phase diagram as 
temperature and pressure are varied22i. At temperatures 
not too low these materials can be regarded as quasi 
one-dimensional. In order to elucidate the relevance 
of bond-charge interactions it is therefore interesting to 
study thermodynamic properties of one-dimensional gen- 
eralized Hubbard models with bond-charge interactions. 
This has been done in ref. |^ for the integrable case t — X 
where t is the single-particle hopping matrix. Here the 
Hamiltonian simplifies considerably since the number of 
doubly-occupied sites is conserved. This allows for an ex- 
act solution for the complete spectruni2242iM which has 
then been used for the calculation of the thermodynamics 
in ref. 21. However, it should be noted that the Hamil- 
tonian at t ^ X has particlc-hole symmetry such that 
the physics is different from Hirsch's scenario. We there- 
fore investigate here the thermodynamics away from the 
point t = X. Since no exact solution is available one 



has to rely on numerical methods. Here we employ the 
transfer matrix DMRG (TMRG) algorithm. 

The TMRG method was originally introduced by 
Nishinc^^ to study two-dimensional classical models. An 
adaption to thermodynamics of quantum systems was 
proposed by Xiang et. al2S*2LS£, who used a Trotter- 
Suzuki decompositionSiffi to map the ID quantum model 
to a 2D classical one. Thereby, various thermodynamic 
properties can be studied with high precision^iiSSiS. 
Similar ideas have recently been applied to stochastic 
modela^^. 

The paper is organised as follows. In Sec. ^ we 
present the model and its known properties. In Sec. IIIII 
the TMRG approach to the thermodynamics of one- 
dimensional quantum systems is explained. Especially we 
emphasize the extensions necessary due to the fermionic 
nature of the problem. Sec . II VI contains our main results. 
We analyze not only thermodynamic quantities like spe- 
cific heat and susceptibilities, but also correlation func- 
tions, expressed by thermal correlation lengths. Sec. [3 
contains a summary and conclusions. 



II. THE MODEL 

The Hamiltonian of the one-dimensional bond-charge 
interaction model reads 

{ij)a {ij)a 

-ft.^(n,| -nil) -/i^(n,-r -hnix) (1) 

i i 

where cl^ and Cia, respectively, are fermionic creation 
and annihilation operators of an up- or down-electron 
(cr =t,i) at site i. In eq. Q the index (ij) denotes 
the summation over adjacent lattice sites, nia- ~ cl^Cia 
the number of particles with spin a at site i and a the 
opposite of spin a. Obviously, the amplitude i > is 
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modified by X > 0, which correlates the hopping process 
to the charge of opposite spins. 

The bond-charge model (Q) exhibits full SU(2) spin and 
U(l) charge symmetry. Under the particle-hole transfor- 
mation 

c,. ^ {~iyc], (2) 

the Hamiltonian transforms (up to constants) into 

H{t,X,n,h) ^Hit-2X,-X,-n,-h) , (3) 

Thus, we can restrict ourselves to the parameter region 
< X/t < 1, which is representative for the whole model. 
< X/t < 0.5 coincides with X/t < such that for weak 
X an repulsive correlated hopping amplitude transforms 
into an attractive one for holes. 0.5 < X/t < 1 conversely 
corresponds to X/t > 1. 

Generically, two universality classes of critical liquid 
phases are found for such two-component systems at tem- 
perature r = 0. In a Tomonaga-Luttinger liquid (TLL) 
phase spin and charge excitations are gapless and all 
correlation functions decay algebraically. In contrast, a 
Luther-Emery liquid (LEL) has a gapped spectrum of 
spin excitations and criticality is only observed for charge 
excitations. 

For X/ 1 <C 1 the bond-charge model (Q) can be treated 
within the framework of bosonisation techniques and 
renormalisation methods^. It was found, that cor- 
responds to an effective Hubbard model with 

C/eff=8XcOs(fcF) , (4) 

and effective hopping amplitude 

toff = (1 - nX)t (5) 

where n is the filling and kp — mi/ 2. Thus, two phases 
appear for X ^ 1 corresponding to the effective Hubbard 
model: a TLL regime (repulsive Hubbard interaction) 
below half filling (n < 1) and a superconducting LEL 
(attractive Hubbard interaction) for n > 1. It is a crucial 
question, if the range of applicability of the bosonisation 
approach holds (qualitatively) for X ^ t. At least we 
know, that the integrable case X/t ~ 1 contradicts with 
the bosonisation results. 

For X/t — 1 the bond-charge model ^ shows an ad- 
ditional SU(2) pseudospin symmetry induced by the so- 
called 77-pair operators'^^ 

i i 

= lT.(^~4^c.1-4lC^l) , (6) 

i 

and is particle- hole invariant due to eq. Q. Apart 
from that, the number of double occupied sites N-^ = 
'^i'l^ifi^il is also conservedS^. 



In that case the model is solvable by mapping the 
Hamiltonian to free spinless fermions22iS^. One can ob- 
serve the complete energy spectrum of Ti. and distin- 
guish three different T — phases. At intermediate 
fillings 0.5 < n < 1.5 some of the ground states are 
of ?7-pairing type, which show ODLRO and hence are 
superconducting^Si2i (see, however, ref. [s^). For small 
fillings n < 0.5 the ground states are identical to those of 
the U — 00 Hubbard model. Due to particle- hole symme- 
try of the model, a corresponding phase occurs for large 
fillings n > 1.5. The model exhibits no spin gap, thus all 
phases fall into the TLL universality class. 

The knowledge of the complete energy spectrum can 
be used to calculate the partition function Z = trc^''^ 
exactlj*2i. This leads to the grand canonical potential 

-13(1) = In (1 + c^'^^^) + - f In (1 + e-'^^"'-''">) dk (7) 

with efc — ~2tcosk and the effective chemical potential 
V M ^ ^ 2cosh/jfe 

The potential (j) can be used to calculate various thermo- 
dynamic properties of the model rigorously^ ^. 

For the non-integrable regime < X/t < 1 detailed 
numerical studies of ground state properties and corre- 
lation functions are available^i^ii^ii^. It was shown that 
the model exhibits a spin gap for < X/t < 0.75 and 
sufficiently large fillings. Conversely, the spin gap closes 
for X/t ^ 0.75 and TLL behaviour is found for all fillings 
n. 

A focus of refs. Il3ll4l was set on various two-point 
correlation functions, such as density, spin, singlet and 
triplet pair correlations. It was found that pair correla- 
tions are dominant for the spin gap regime < AT ^ 0.75. 
X/t sa 0.5 was identified as a particular point where the 
spin gap and the dominance of superconducting correla- 
tions are maximal. 

The present paper expands the T = studies to fi- 
nite temperatures T > 0. Apart from the investigation 
of susceptibilities and the specific heat, we particularly 
focus on thermal correlation functions, especially the pos- 
sibility of dominant superconducting pair correlations at 
T > 0. For the numerical analysis we use the transfer 
matrix DMRG (TMRG) method, which has been shown 
to be highly precise for fermionic modcls'^^'^'^ without 
suffering from minus-sign problems etc. There it has 
also been demonstrated, that the asymptotic behaviour 
of various correlation functions in terms of thermal cor- 
relation lengths are accessible by the TMRG algorithm. 

III. THE TMRG METHOD 

The transfer matrix DMRG (TMRG) method, a vari- 
ant of White's DMRG algorithm^^'^o^ has originally been 
introduced by Nishino to study two-dimensional classical 
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FIG. 1: Suzuki- Trotter decomposition of the partition func- 
tion Z, leading to a classical two-dimensional lattice. The 
dots (•) represent the lattice sites of classical spins, which 
interact by the shades plaquettes. In (a) the traditional map- 
ping is plotted, cf. ref . 1301. An alternative lattice is depicted in 
(b), which has additional auxiliary lattice points, cf. ref. 
The quantum transfer matrix Tm is constructed column-wise. 



systemsS^. Xiang et. al applied the TMRG to the quan- 
tum transfer matrix (QTM) to study the thermodynam- 
ics of one-dimensional quantum systemsSSi2L2^. Since 
then it has been successfully used to study various spin 
and fermion systems^ ^■12*^. 



A. Transfer matrix approach 

The key idea of the TMRG is a Suzuki- Trotter 
decompositio n^^i^" of the partition function 



Z = tr c 



= trie 



-«°e-^«")*' + 0(e2) (9) 



where e = (3/M is fixed and has to be chosen suf- 
ficiently small. The Hamiltonian consists of next- 
neighbour interactions Ti — and is spht up into 
7io/o = So/c^«,j+i' where the sum is taken over odd 
and even bonds, respectively. The decomposition leads 
to a two-dimensional classical model with a checkerboard 
structure (Fig.n(a)) and plaquette interactions r given 
by 



One finds that the thermodynamics in the thermody- 
namic limit L — > c» is determined by the leading eigen- 

I L / R\ 

values Aq and eigenvectors \Aa ) (a = 0, 1, . . . ) of Tm- 
Note, that left and right eigenvectors have to be distin- 
guished because Tm is in general not symmetric. The 
free energy explicitly reads 



/ = -|lnAo 



(11) 



where the leading eigenvalue Aq is real and unique due 
to Frobenius' theorem. Thus, we always observe a gap in 
T]\j , which only may close and lead to a phase transition 
for M — + oo (T ^ 0) . The expectation value of any local 
operator Oj is given by 



{A^\Tm{Oj)\A^) 
Ao 



(12) 



where in TM{Oj) one plaquette of the QTM Tm has been 
modified by the operator Oj. 

In the TMRG algorithm, the matrices Tm and TM{0-j) 
are sequentially enlarged using a DMRG-like algorithm. 
An asymmetric density matrix p = tr' |A^) (AqI is used 
to target the leading eigenvector. We do not discuss the 
details, since the TMRG we use basically follows Ref.l2^ 
[2I. As spin and particle number are conserved by 7i, 
Tm has a block structure that drastically reduces the 
computational effort. Alternatively, we use the number 
of particles with spin up and down TV^. = riicr to build 
the QTM blocks. To be more precise, N„ are conserved 
for a row-transfer matrix. Contrary, lS.N„ — X]i(~l)'"'*'^ 
are the corresponding good quantum numbers for the 
(column) quantum transfer matrix Tm<^ Hence, we label 
the QTM block by the tupel (AiV^, AiV^). The leading 
eigenvalue Aq is located in the (0,0) block. This follows 
directly from the high temperature limit T — s- 00, where 



-phi 



is unity and therefore 



Jiri 



hr2 ^hri 



(13) 



Thus, only matrix elements of the (0, 0) block of Tm 
are non-vanishing. Consequently Aq is also found in the 
(0,0) block. This property persists for finite tempera- 
tures T > 0, because no level crossing is expected, which 
would imply a phase transition. 



{hi 
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\rir2) 



(10) 



The quantum direction L is thereby supplemented by a 
virtual Trotter dimension M . The plaquettes are column- 
wise combined to the quantum transfer matrix (QTM) 
Tm, cf. Fig. nja). Alternatively, the decomposition of 
the partition can be constructed^ as shown in Fig.^b), 
which is fully translational invariant in quantum direc- 
tion, in contrast to (a) . The advantage of (b) is discussed 
later, we first focus on variant (a). 



B. Correlation lengths 

One focus of the present work are correlation functions, 
such as spin, density and superconducting pair correla- 
tions. The general form of a two-point correlation func- 
tion (OjOj) for a local operator Oj at finite temperatures 
reads 



{oIo,)-{oI){o,) = J2m^ 



(14) 
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All correlations decay exponentially because no phase 
transition can occur for T > 0. The correlation lengths 
(CL) are directly accessible by the leading eigenvalues 



.-1 1 1 Ao 



(15) 



The wavevectors ka are given by the complex argument 



1 



Ao 

Aa 



+ riTT {n — 0, 1) 



(16) 



As the classical lattice in Fig. ^a) is translational in- 
variant in quantum direction, but with a doubled unit 
cell, ka is only determined modulo tt. This problem can 
be solved by using the QTM presented in Fig. ^b) , cf. 
ref.lH 

The crucial question which eigenvalue Aq, refers to 
which correlation function is controlled by the coefficients 
Ma, that are given bySS 



Ma = 



{A^\TMiOl)\A^){A^\TM{Oj)\A§) 

AoAa 



(17) 



Eq. ifTTI) represents a selection rule to assign A^ to the 
correct type of correlation length. A first criterion are the 
quantum numbers (AiV|,A7V|) of jAa^'^): Ma IS non- 
vanishing only, if (A^l and TuiOj) \Aq) have the same 
quantum numbers. To be more concrete, we discuss now 
the relevant correlation functions. 

Longitudinal spin {Oj — nj^ — nji) and charge {Oj — 
riji correlations are determined by eigenvalues A^ 

of the (0, 0) block of Tm, since Oj is diagonal and TM{Oj) 
does not change the quantum numbers. They have to be 
distinguished by calculating Ma explicitely. 

The situation changes for the superconducting singlet 
and triplet pair correlations, where Oj in eq. (|14() is given 

by 



ptl 
j 



V2 

1 



(18) 

(19) 
(20) 



The pair correlation operators change the quantum num- 
bers. Therefore, eigenvalues A^ of the (±1,±1) block 
contribute to the correlation lengths of Pj and P*^, 
whereas those of the (±2, 0) and (0, ±2) block belong to 
Pj^ and Pj^, respectively. Without magnetic field, the 
triplet correlation functions Pp are equal. Thus one can 
unambiguously distinguish singlet and triplet pair corre- 
lation lengths by computing the eigenvalues Aq of the 
(1, 1) and (2, 0) block only. 

Note, that the present work uses both types of QTMs 
of Fig. ^ Even if the second variant improves the calcu- 
lation of ka, the checkerboard style QTM (Fig. QJa)) is 
more precise, especially for the correlation lengths. The 
reasons for that are technical ones within the TMRG al- 
gorithm. 



C. Jordan- Wigner Transformation 

A sophisticated point in the TMRG algorithm are 
the fermion statistics which have not been considered in 
the Suzuki- Trotter mapping in eq. Hence we use 

a Jordan- Wigner transformation (JWT)"'^ to map the 
fermion system onto a spin model. We define two unitary 
operators 



K 



exp 



K 



exp I ittT^ StS„ 



+ a- 



(21) 



(22) 



where 5". (5^^) are spin-1/2 lowering (raising) matrices 
of two "spin types" a =t,i at site i. These are related 
to the fermion operators by 

c,,=K,,S~„ and c]„ = S%K,„. (23) 

These operators fulfill the Fermi statistics. Inserting the 
JWT into the Hamiltonian the hopping terms cjcj 
modify, whereas the diagonal charge operators jii^ trans- 
form canonically: 



"1+ 



(24) 



Note that periodic boundary conditions may transform to 
twisted ones, but this is irrelevant for our studies where 
the thermodynamic limit is performed exactly. We point 
out that without the JWT the TMRG algorithm com- 
putes wrong results in particular for the singlet correla- 
tion lengths of the (±1, ±1) block, although many other 
quantities are not affected. 



IV. RESULTS 

We present numerical TMRG results for the parameter 
region < Xjt < 1, which is representative for all X/t, 
cf. Sec.|nJ In all computations we have set t = 1 and tem- 
peratures are measured in units of t. As we use a grand 
canonical description, the chemical potential fj, controls 
the density n of particles. Note, that the density is there- 
fore temperature dependent, n = n{T,fi). If not stated 
differently, we use e = 0.05 and retain = 70 — 100 
states in the TMRG algorithm. This is already suffi- 
cient for highly accurate results. We typically compute 
M « 1000 TMRG iterations, which correspond to a min- 
imal reachable temperature of T = 0.02. 



A. Thermodynamics 

In this section we analyze various thermodynamic 
quantities: the grand canonical potential 0, local ex- 
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pectation values (e.g. local density n), charge and spin 
susceptibilities Xc and Xs and the specific heat c^. 

The grand canonical potential (^(T, ji) is directly ac- 
cessible through the largest eigenvalue of the QTM, cf. 
Sec. mil The charge susceptibility 



dn I 



(25) 



is obtained by computing the density n = (n^i + riii) of 
electrons for two different chemical potentials fi — Afi and 
fj, + A/i. Then, we use numerical derivatives 



Xc(r,pi) 



n{T, II + A/i) — n(T, /i — A/i) 
2A/i 



(26) 



to approximate Xc where typically A/i ^10 ^ is chosen. 
Similarly, the spin susceptibility 



h=0 



dm I 



(27) 



is obtained by the numerical derivative of the magneti- 
sation m = (riii — nil), varying the magnetic field h. 

The natural form of the specific heat in a grand canon- 
ical formulation is given by 



T 



ds I 



T 



I p 



(28) 



where s denotes the entropy per site. But usually one 
finds data for the specific heat at constant density n, 
which is related to by the thermodynamic relation 



T 



ds I 



T 



dn , dfi. 



(29) 



We only determine by calculate numerically the second 
derivative of the grand canonical potential (p. 



1. Exactly Solvable Case X = 1 

First, we show data for the exactly solvable case X = 1, 
basically to check the precision of the TMRG algorithm. 
In Fig. 12 the grand canonical potential 4> is plotted for 
the chemical potentials /i = —0.5,0,0.5. The tempera- 
ture dependence of the density n(T,fi) is plotted in the 
upper inset. The particle-hole symmetry of the model 
is confirmed by the symmetry of n(T, ±0.5) as well as 
by the potentials (/)(T, ±0.5), which are only shifted by 
A/t = 1. Our data are compared to exact results, calcu- 
lated by eq. Q, cf. the lower inset of Fig. [3 The absolute 
error is less than 10^^ for temperatures T > 0.1 and less 
than 10~^ for 0.02 < T < 0.1. Not explicitly shown here, 
the errors of the density n as well as other expectation 
values (e.g. the magnetisation) are approximately of the 
same order. 

The origin of a numerical error of at least 10^^ can 
be understood by the approximation done in the Suzuki- 
Trotter decomposition in eq. (|5J|. The choice of e (in our 
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FIG. 2: (left) Grand canonical potential (p as function of tem- 
perature for X — 1 and chemical potentials /i = —0.5, 0, 0.5. 
The absolute deviation from the exact value in dependence of 
the temperature is plotted in the left inset for jj. — 0. The 
corresponding TMRG data for the density n = n(T, n) is de- 
picted in the right inset. For comparison exact data using 
eq. lO are shown by circles (o). 
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FIG. 3: Density n(/i, T) as a function of temperature and 
chemical potential for X = 1. For T — > a discontinuity at 
/i = is found, where the density jumps from n = 0.5 to 
n = 1.5. 



case e = 0.05) causes an error of the order 0{e^) in the 
partition function l^. 

It is interesting to investigate also the density n{T, fi) 
as a function of /i and T (Fig. OJ. For T ^ one ob- 
serves the formation of a plateau for /i = 0. This can be 
explained by the ground states of the 77-pair phase, which 
occurs for 0.5 < n < 1.5. Here, all ground states with 
different numbers of pairs are degenerate. Thus the 
chemical potential does not change by adding a particle, 
because a new bound pair is build without changing the 
energy. Thus, the thermodynamics of the 77-pair phase 
can be only realized for /i = 0. Any /i 7^ induces n < 0.5 
or n > 1.5 and falls into the U — 00 phase, cf. Sec. ITU 
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FIG. 4: Spin and charge susceptibility as a function of tem- 
perature for X = 1 and ^ — 0,0.5, 1.0. The upper diagram 
plots Xa ■ T over T in order to demonstrate that Xs diverges 
like Xs '~ 1/T. The insets depict the relative deviation from 
exact results for ^ = 0. 
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FIG. 5; Spin susceptibility as a function of temperature for 
various amplitudes X/t = 0.9... 0.1 and fixed chemical po- 
tential \i = 0.6. The upper figure shows Xs for X — 0.9 . . . 0.5 
the lower one for X = 0.1 . . . 0.5, respectively. The inset de- 
picts the density n{X,fi = 0.6, T) for various parameters X. 



Fig. 2] depicts the spin and charge susceptibilities for 
X = 1 and /i = 0,0.5,1.0. The comparison with exact 
data results in a relative error less than 10^^ for temper- 
atures T > 0.1. Not surprisingly, it is larger than that 
for 0, because numerical errors are propagating through 
the numerical derivatives. 

From conformal field theory (CFT) the low tempera- 
ture limit T ^ of and Xc are respectively given by 

2 2K 

Xs{T~^0) = and Xc{T ^ 0) = — , (30) 

TTVs TTVc 

where Vs and Vc denote the spin and charge velocities and 
Kc the LL parameter. This can be used to interpret the 
results of Fig. 0] The ground states of X = 1 for n < 
0.5 (and after particle-hole transformation for n > 1.5) 
corresponds to the U = oo Hubbard model, cf. Sec. 
For the U = oo Hubbard model it is known, that 

Vs = and Vc — 2t|sin(7rn)| . (31) 

The spin excitations exhibit no dispersion due to com- 
plete degeneracj*2i. In Fig. 0] we have plotted Xs • ^ to 
demonstrate, that the spin susceptibility Xs is diverg- 
ing for all /i and fillings n in the limit T ^ since 
limT^o Txs is finite. The charge susceptibility Xc is di- 
vergent for /i = 0. 

2. Non-Integrable Case < X < 1 

We now consider the non-integrable case < X < 1. 
From analytical-'-' and numericaU^ii^*i^ approaches it is 



known, that a crossover from a TLL phase to a super- 
conducting LEL is observed. For X <C 1 this is evident, 
because the bond-charge model is an effective Hubbard 
model (cf. Sec. where the crossover is expected to 
appear at UoS = 0, corresponding to half- filling n — 1 
(see Q). 

In our thermodynamic study we can identify the exis- 
tence of a spin gap if the susceptibility Xs vanishes for 
T 0. In Fig.|Slxs is depicted for fixed chemical poten- 
tial = 0.6 and various hopping amplitudes < X < 1. 
For X > 0.65 we can clearly verify TLL behaviour, be- 
cause the susceptibility does not vanish for T — > 0. At 
X « 0.65 we observe a crossover to a LEL phase with 
Xs{T ^ 0) = 0. Decreasing X further, the gap seems to 
close again for X < 0.3. 

Additional the inset of Fig. |S1 shows the density 
ri(r, /i = 0.6) for the respective parameters X. Note, 
that fi — 0.6 induces a regime of more than half filling 
n > 1 for each X. Therefore the spin gap should per- 
sist especially for AT <C 1 , which we can not show by our 
TMRG results. But from the attractive Hubbard model 
it is known, that the spin gap is exponentially small for 
t/ < 1. Since from ^ we have C/off < 1 for A < 1, the 
gap in the bond-charge model becomes extremely small, 
especially at half filling. The decay of Xs at A ^ 1 is 
therefore expected at such low temperatures that are not 
accessible by TMRG. To verify that aU A < 0.5 exhibit 
a spin gap we increased the filling in Fig. 1^1 Here, the ef- 
fective Coulomb potential J7cff is larger and the spin gap 
appears in our data. 

In order to give a more detailed overview of the ther- 
modynamics of the bond-charge model for < A < 1, 
we focus on three particular points A = 0.1, 0.5, 0.9. 
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Temperature T 

FIG. 7: Spin and charge susceptibilities Xs and Xc as a func- 
tion of temperature for X/t — 0.1 and various chemical poten- 
tials ^. The inset plots the density n{fj,,T). The correspond- 
ing T = values for the Hubbard model with effective Uctf 
and teff are shown by symbols, (o) : jj. — —1.5, (□) : /i = —1, 
(O) : ^ = and (A) : = 1.5. 



Fig. Ul plots the charge and spin susceptibility a± X = 
0.1 for various chemical potentials /x. Additionally, in 
the inset, the density n{ii,T) is depicted. As X is small 
compared to the bandwidth t, the model iQJ should co- 
incide an the effective Hubbard model. We check the 
asymptotics of the susceptibilities Xs and Xc for T — > 
using eq. (|Sn|l . These can be calculated exactly for the 
integrable Hubbard model^Sii^. For a weak (repulsive) 



FIG. 8: Specific heat as a function of temperature for 
X = 0.1 and chemical potentials /i < (upper figure) and 
/i > (lower figure). 



Coulomb potential J7, one finds 



Vc/s = Vp[\± 



. (32) 



U 



Ant sm{kp) 



(33) 



with VF = 2tsin(fcF) and kp = Using eqs. (ESJ 
and (|30|l we have calculated Xs and Xc for T ^ (see 
Fig- El- The density no = n{fi,T 0) which was used 
to determine kp, UeS and tes is denoted in the legend 
of Fig. 13 We observe perfect agreement with our data, 
which again supports the correspondence of (effective) 
Hubbard and bond-charge model. 

Even more evidence of Hubbard like thermodynamics 
is given by the shape of the specific heat . As shown for 

in Fig. |S1 two characteristic features can be observed 
for large and small fillings, a shoulder at low tempera- 
tures and a peak at a slightly larger temperature. As in 
the repulsive Hubbard model these features can be re- 
lated to spin and charge excitations, respectively. They 
merge close to half filling where we have effectively a free 
fermion system for X ^ 1 . Moreover the regime less than 
half filling n < 1 closely corresponds to n > 1. Again, 
this can be explained by the exponentially small gap, 
which does not affect the physics at those temperature 
scales we can observe by TMRG. 

The effect of the spin gap manifests itself at larger 
interactions X. Fig. |51 depicts charge and spin suscep- 
tibilites for X = 0.5. In contrast to X = 0.1, the spin 
susceptibility Xs clearly affirms a gap for n > I since 
Xs{T ^ 0) — > 0. Therefore the phase transition oc- 
curs at half filling, which is predicted by the bosonisation 
results^. A significantly larger spin gap than at X = 0.1 
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FIG. 9: Spin and charge susceptibilities Xs and Xc as a func- 
tion of temperature are depicted for X — 0.5 and various 
chemical potentials /j.. The inset plots the density n{^,T). 
The corresponding T — values for Xs and Xc calculated from 
the Hubbard model with effective UcS and tcs are shown by 
symbols, (o) : fj, ^ —1, (□) : ^ = 0, (O) : ^ = 0.4 and 
(A) : /i = 0.5. 



tl=-2 

H=-1.5 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Temperature T 

FIG. 10: Specific heat Cfi as a function of temperature for 
X = 0.5 and various chemical potentials fi < (upper figure) 
and n > (lower figure). 



FIG. 11: Spin and charge susceptibilities Xs and Xc are de- 
picted for X = 0.9 and various chemical potentials /i. The 
inset shows the density n{fi,T). 



X ~ 0.5. Since Ucs is not small here, we used the Bethe 
Ansatz (instead of eq. H32|l and H33() ) to calculate the 
T ^ values of Xs and Xc for the (effective) Hubbard 
model, which are also depicted in Fig.lHl Obviously, these 
do not fit to our TMRG data. 

For large X w 1 the spin gapped phase disappears, 
which also contradicts the results of bosonisation. As an 
example we show TMRG data for X — 0.9. The sus- 
ceptibilities and the density, which are shown in Fig. 1111 
qualitatively coincide with the X ^ 1 case. At /i « 
we observe a jump in the density, indicating a phase 
comparable to the ry-pair phase. The spin susceptibil- 
ity XsiT — > 0) is diverging for all fillings n, which yields 
Vs 0. Fig. El additionally plots the specific heat 
for various chemical potentials /i. In fact the model is 
observed to be nearly particle-hole symmetric {fi — /z). 
Note, that for = the specific heat exhibits an ad- 
ditional low energy peak. Spin excitations can not con- 
tribute to that peak due to complete degenerac y. For 
X = 1 such peak has also been observed in ref. |2ll It 
is associated with the melting of pairs. Remarkably, the 
physics at A" = 0.9 is very similar to the highly symmet- 
ric A = 1 case. Thus, A = 1 is not a singular point, but 
the characteristics persist for a certain neighbourhood 
around A = 1. 



is also verified by the specific heat shown in Fig.^l The 
linear behaviour of Cfj, for fj- ^ 1 occurs only at very low 
temperatures and can not be observed on the tempera- 
ture scale shown. Here the exponential corrections com- 
ing from the finite spin gap dominate. 

Even though the numerical data for the thermody- 
namics suggest Hubbard like behaviour on a qualitative 
level, a more detailed quantitative comparison fails for 



3. Universal Crossing Points 

In the thermodynamic data shown above we find vari- 
ous phenomena of (nearly) universal crossing points. As 
a first example, the specific heat in Fig. |S1 exhibits even 
two crossing points. But also in the density in Fig. [3 or 
Fig. 13 a crossing point is apparent at n w 1. 

In the case of the Hubbard model, it is well known^ 
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1.0 1.5 2.0 

Temperature T 



FIG. 12: Specific lieat as a function of temperature for 
X = 0.9 and various cliemical potentials /i < (upper figure) 
and /i > (lower figure). 



that such crossing points appear. In refs. I45ll46l general 
thermodynamic arguments for the existence of a sharp 
crossing point in the specific heat are given. 

Sharp universal crossing points {nc,nc) in the density 
n =^ n{T, fj.) can be generally explained in a similar way. 
The integral 



On I 

5^1/^ = 



°° ds. ^ 



s(oo, T) 



s(-cx),T) = (34) 



vanishes: for /i — )■ ±oo the band is either completely 
empty or filled, thus s ^ 0. Then has to change 

its sign. Due to continuity arguments we conclude the 
existence of fic, such that 



dn 



n{T, fic) — const. 



(35) 



In the case of the bond-charge model at weak coupling 
X <C 1, the crossing point of the density appears at half 
filling 71c = 1 , where the model effectively behaves as free 
fermions {Ucs — 0). 



B. Correlation Lengths 

A great advantage of the TMRG algorithm is that 
not only quantities related to the free energy Hll|l are 
accessible, but also, through the next-leading eigenval- 
ues of the transfer matrix, thermal correlation lengths 
(|15|l . As explained in Sec. IIII 51 the corresponding cor- 
relation functions are identified by their quantum num- 
bers (AN^i^ANi). It is even possible to distinguish 



the contributions of different wavevectors ka = akp 
[a ~ 0, 2,4, . . .). This will be done in the following to 
determine the dominant correlations in the different pa- 
rameter regimes. First we discuss the low-temperature 
behaviour, which can be compared to predictions of con- 
formal field theory. 



1. Low-temperature behaviour 

The correlation functions of a TLL at temperature 
T — Q show universal behaviour. According to confor- 
mal field theory (CFT), density-density (d-d), spin-spin 
(s-s), singlet pair (sp) and triplet pair (tp) correlation 
functions readii^ 



(n^no) 



{p^'n) 



Aor-^ + ^ir-(i+^'^) cos(2fci.r) 
Bqv-^ + Bir-''^+^-='^ coa{2kFr) 

+Cir"^^'^+"^)cos(2fcFr-) , 



Dnr 



-(i+ifc) 



+D,r-(''^+Tt+^) cos(2fc. 



(36) 
(37) 

(38) 

(39) 



(up to logarithmic corrections). Thus, all critical expo- 
nents are controlled by a dimensionless parameter Kc- 

In contrast, in a LEL the spin excitations are gapped. 
Thus s-s and tp correlations decay exponentially and do 
not show universal behaviour. According to CFT, the 
asymptotic behaviour of the d-d and sp correlations is 
given by^^ 



{nrUo) ^ Aor ^ + Air cos{2kFr) 

+A2r-^'^- cosiAkpr) , (40) 

+Cir-'-^-'+Thhos{2kFr) . (41) 

In the present work we are discussing correlation func- 
tions at finite temperatures T > 0. Here, the model is 
always non-critical and all correlation functions decay ex- 
ponentially. The asymptotics are described in terms of 
thermal correlations lengths (CL), cf. eq. I|14|l . which are 
temperature dependent. 

The asymptotic behaviour of the CL for small temper- 
atures T ^ can stiU be obtained by CFT^^. For a TLL 
all CLs diverge 



e(T ^ 0) 



1 
T 



(42) 



The scaling dimensions Xs and Xc can be calculated ex- 
plicitly. One finds, that the non-oscillating d-d (s-s) CL 
is given by 



(0) 

7dd 



2tt' 



(0) 



Vs_ 

2tt 



(43) 
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The 2kp parts of both read 



(2fcf) 
'dd/ss 



whereas the 4k p part of the d-d CL yields 

Ukp) Vc 

7dd - 



(44) 



(45) 



The non-oscihating parts of sp and tp CL are given by 

Vc 



(0) 
'^sp/cp 



the 2kp part of the sp CL by 



/sp 



if. 



and the 2kF part of the tp CL by 

(2kF) Vc 

7tp - 



^l^+Kc + 2^ 



(46) 



(47) 



(48) 



In the LEL the tp and s-s CLs do not diverge for T 
due to the spin gap. But the asymptotic behaviour 
of the d-d and sp CLs is still predictable by CFT and 
corresponds to eq. H43|l - H47|) in the limit Us — > oo. Thus, 
we have 



JO) 



^dd ~ 27r' '^'^'^ 



(2fcF) 



and 7^^) 



for the non-oscillating, 2kp and 4kp d-d CLs and 

Vr 



Tsp 



(0) = !^ and 7if-) 



K, 



(49) 
(50) 



for the non-oscillating and 2kF sp CL, respectively. Note, 
that the low temperature asymptotics of the correlation 
lengths depend not only on the universal exponent Kc, 
but additionally on the charge and spin velocities Vc and 



2. Comparison of correlation lengths 



It was shown in Sec. lIIII that the TMRG algorithm pro- 
vides a very effective facility to compute thermal corre- 
lation lengths numerically. First, we check the precision 
of our TMRG data. For that purpose we choose X — 0, 
which describes a system of free fermions with spin. Here, 
all correlation lengths can be computed exactly. Due to 
Wick's theorem, spin, charge and pair CLs are identical 
given hy^ 



= 2 arsinh(7rT/2) 



(51) 



Fig. El compares the CLs computed by the TMRG 
program to exact data. Down to a temperature T sa 0.2 



10" 



10 ■ 



10 ' 



10" 



10 



density-density 
triplet pair 
singlet pair 




2 3 
Temperature T 



FIG. 13: Precision check of various correlation lengths com- 
puted by the TMRG algorithm for free fermions. The relative 
error is plotted as a function of temperature T, showing that 
the GLs keep reliable for T > 0.2. 



the relative error is shown to be less then 10%. As already 
mentioned in Sec. Illll it is important to perform a JWT 
before applying the TMRG to the bond-charge model. 
Otherwise, particularly the singlet pair correlations are 
not correctly reproduced by the algorithm. 

We now discuss in detail s-s, d-d, sp and tp correlations 
in the bond-charge model. We focus on the parameter 
point X = 0.5, where the spin gap is comparatively large 
for n > 1. For T = 0, the system is then a LEL with dom- 
inating superconducting sp correlations due to Kc > li^ 
From eqs. ijl^ and it is expected, that sp correla- 
tions should dominate even for finite low temperatures 
T > 0. At zero temperature, close to half- filling a tran- 
sition to a TLL takes place which exists for all densities 
n < 1. 

Fig. 1141 shows the (leading) thermal CLs for three dif- 
ferent chemical potentials (/i — —1,0.4,0.6). The choice 
of /i covers the TLL (/z = —1) and the LEL phase 
(/i — 0.6) as well as a point close to the phase transition 
(/i = 0.4). In these figures we have plotted ^ ■ T vs. T 
which allows to determine the low temperature asymp- 
totics. As ^ = 7/r for T ^ (cf. eq. the curves 
should become linear and intersect the ordinate at 7. Ac- 
cording to CFT (see Sec. lIV iHTl . the factor 7 depends on 
the spin and charge velocities and the critical exponent 
Kc. 

For /i = —1 (no ~ 0.38) in Fig. E| (a) the system is 
in the TLL regime. The non-oscillating d-d CL domi- 
nates for all temperatures. The leading s-s CL for T ^ 
shows incommensurable oscillations (fc = 2kp), whereas 
the k = part is strongly suppressed and not shown in 
the figure. Additionally, the k — 4kp d-d correlation 
lengths dominate the k = 2kp d-d and all s-s CLs. This 
scenario can be easily understood by the fact, that the 
spin velocity is very small (u^ <C Wc), which can be ver- 
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spin-spin {k=2kp) 
density-density (/:=0) 
density-density (/:=4/:^) 
triplet pair (A:=0) 
singlet pair (k=2kp) 



Temperature T 



spin-spin (k=2k^) 

density-density (k=0) 

density-density (k=4kj^) 

triplet pair (A:^=0) 

singlet pair (^=2^^) 




0.4 
^ 0.2 



0.0 
0.4 

^ 0.2 



- (a) 



0.0 
0.4 



^ 0.2 



density-density (k=0) 
singlet pair (^=0) 
singlet pair {k=2k^) 



1 1 1 1 1 — 

- (b) 


T ' 1 ' 1 ' 1 ' 1 ' 1 1 1 1 

/; — 

1 . 1 . 1 . 1 . 1 . 1 . 1 . 


1 . 1 . 1 . 

- (c) 


1 . 1 . 1 . 1 . 1 . 1 . 
'/; ^^^^-^-'^ 



0.0 0.2 0.4 0.6 



0.8 1.0 1.2 
Temperature T 



1.6 



1.8 2.0 



FIG. 15: Crossing phenomena of density-density (fc = 0) 
and singlet pair (k = 0,2kF) correlation lengths for various 
chemical potentials (a) jj, — 0.6 (b) n — 0.8 and (c) fil.O. 
The diagrams plot the particular correlation length ^ ■ T as 
a function of temperature T. The corresponding densities 
no = n{T — » 0) are given by (a) no ~ 1.21 (b) no ~ 1.92 (c) 
no « 1.99. 



Temperature T 



(c) 






spin-spin {k=2k^ 






density-density (/;^-=0) 






density-density {k=Ak^ 






triplet pair (^^=0) 






singlet pair (/:^=0) 






singlet pair (k=2kj 









0.0 



1.0 2.0 

Temperature T 



3.0 



FIG. 14: Plot of the leading d-d, s-s, sp and tp correlation 
lengths ^ for X — 0.5 and (a) fi — —1, (b) /i = 0.4 and (c) 
fi = 0.6. The corresponding densities no = n{T 0) read (a) 
no ~ 0.38, (b) no ~ 0.88 and (c) no = 1.19. The diagrams 
show the respective correlation ^ ■ T as a function of T. As 
an example, the inset of (a) depicts the wavevectors k of the 
oscillating CLs. The the circles (o) correspond to the T — > 
limit of fc = 2fcF and k — Akp with kp ~ nn/2. 



which agree perfectly with the TMRG data. 

The case = 0.4 (no w 0.88) shown in Fig. [H (b) 
also falls into the TLL regime, but is situated close 
to the phase transition to the LEL. The CLs generally 
get smaller, because the charge velocity Vc is decreased 
(Fig. I^J. Note, that the crossing phenomena of CLs be- 
come very rich close to the transition. 

A crossover of sp and d-d correlations at finite tem- 
perature Tc « 0.5 is observed in Fig. (c) for /i — 0.6 
(no « 1.19). And in contrast to the latter cases, the lead- 
ing sp correlation length is shown to be commensurable. 
This agrees with the predictions of eq. 149() and H50() . if 
i^c > 1 is assumed. 

The crossing temperature Tc increases for higher fill- 
ings. This is demonstrated by Fig. ^1 which shows 
only the leading d-d and sp correlation lengths for /i = 
0.6, 0.8, 1.0. Note, that the CLs in Fig.UKc) and Fig.lTsl 
do not show asymptotic ^ behaviour for the achievable 
temperature region T > 0.2. 



CONCLUSIONS 



ificd by Fig. ^ The ratio Vc/vs ^ 1 enters eq. (Hit , 
thus these correlations are suppressed for low tempera- 
tures. For the same reason, the tp correlation length is 
crossed over by the leading incommensurable sp correla- 
tion length at low temperatures, cf. eq. H46|l . The inset 
of Fig. El (a) depicts the wavevectors kp. We have also 
plotted the T — > values of 2kp and Akp with kp = 



In this paper we have studied in detail the thermody- 
namics of the bond-charge model by means of numerical 
TMRG calculations. We focused on the parameter re- 
gion < X/t < 1, which is representative for the whole 
model due to particle-hole symmetry. 

Our computations show data for the grand canonical 
potential, charge and spin susceptibilities, the specific 
heat, density and thermal correlations lengths. The de- 
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tailed study of small {X/t = 0.1), intermediate {X = 0.5) 
and large {X = 0.9) bond-charge coupling was performed 
to determine the underlying physics of the model. 

In accordance with previous works, numerical studies 
of the spin and charge susceptibilies exhibit two T = 
phases of the model. For less than half-filling we prin- 
cipally find Tomonaga-Luttinger liquid behaviour, where 
spin and charge excitations are gapless. For more than 
half-filling a Luther-Emery liquid phase with spin gap is 
found, if X is not too large. 

We have shown that the thermodynamics for small 
X <^ 1 essentially coincide with that of the Hubbard with 
an effective Coulomb potential J/eff and hopping ampli- 
tude teS, as predicted by bosonisation. For intermediate 
coupling (X = 0.5) the correspondence holds qualita- 
tively, but not quantitatively. For large X « 1, the spin 
gap disappears for all fillings, which even qualitatively 
contradicts the bosonisation results. 

Additionally, numerical investigations of the specific 
heat and density exhibit interesting phenomena of nearly 
universal crossing points. These can be understood from 
quite general thermodynamic considerations. 

Finally we focused on thermal correlation functions, 
which we characterized by thermal correlation lengths. 



In the Luther-Emery liquid phase, we observed interest- 
ing crossing phenomena of superconducting singlet-pair 
and density-density correlations at finite temperatures 
Tc > 0. Even though true superconductivity is not possi- 
ble in the one-dimensional bond-charge model, our stud- 
ies support a strong tendency towards superconductiv- 
ity. This indicates, that the bond-charge interaction may 
be relevant for superconducting phenomena in higher di- 
mensions. 

As we concentrated on the physics of correlated hop- 
ping only, additional interactions like on-site Coulomb re- 
pulsion U or pair hopping processes Y were not discussed 
in this paper. Work is in progress in these matters. 
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